function sse = FitLorentzian(p, xdata,ydata )
        FittedCurve = lorentzian(p,xdata);      
        ErrorVector = FittedCurve - ydata;
        sse = sum(ErrorVector .^ 2);
end